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1 Introduction 



We consider exactly solvable models of statistical mechanics in one space and one time dimension. 
The Quantum Inverse Scattering Method and Algebraic Bethe Ansatz are effective methods for a 
description of the spectrum of these models. Our aim is the evaluation of correlation functions of 
exactly solvable models. Our approach is based on the determinant representation for correlation 
functions. It consists of a few steps: first the correlation function is represented as a determinant 
of a Fredholm integral operator, second the Fredholm integral operator is described by a classical 
completely integrable equation, third the classical completely integrable equation is solved by means 
of the Riemann-Hilbert problem. This permits us to evaluate the long distance and large time 
asymptotics of the correlation function. The method is described in 0. The most interesting corre- 
lation functions are time dependent correlation functions. The determinant representation for time 
and temperature dependent correlation functions of quantum nonlinear Schrodinger equation was 
obtained in . In this paper we describe the correlation function by means of completely integrable 
integro-differential equation. In the forthcoming publication we shall formulate the Riemann-Hilbert 
problem for this equation and evaluate long distance asymptotic. 

The quantum nonlinear Schrodinger equation can be described in terms of canonical Bose fields 
ip{x,t),'ilj^ {x,t) {x G R) obeying 

[V'(x,t),V^(y,i)] =(5(x-y). (1.1) 
The Hamiltonian and momentum of the model are 



H 



dx (^dxijj'^ {x)dx'4^{x) + ciIj\x)iIj\x)iI){x)iIj{x) — hip\x)'tp{x)^ , (1-2) 



P = -i J dxi;^{x)dxilj{x). (1.3) 
Here < c < oo is the coupling constant and /i > is the chemical potential. The spectrum of the 



model was first described by E. H. Lieb and W. Liniger [13|, |14|. The Lax representation for the 
corresponding classical equation of motion 

i^^i;=[i;,H] = -^^|; + 2c^|J^^-h^P, (1.4) 



was found by V. E. Zakharov and A. B. Shabat [17|. The Quantum Inverse Scattering Method 
for the model was formulated by L. D. Faddeev and E. K. Sklyanin [Q]. The quantum nonlinear 
Schrodinger equation is equivalent to the Bose gas with delta-function interaction. In the sector 
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with particles the Hamiltonian of Bose gas is given by 

N q2 

^^ = -E^ + 2c E S{zk-Zj)-Nh. (1.5) 

j=l '^^j i-<j<k<N 

In this paper we shall consider the thermodynamic of the model. The partition function and the 
free energy of the model are defined by 

Z = tre^^=e~^. (1.6) 
The free energy F can be expressed in terms of Yang- Yang equation [16] 

oo 

eiX) = A^-.-i:/-,_^|_^,n(l + .-r),,, (1.7) 

— oo 

OO 

F = J In (^1 + e-'-^^ dfi. (1.8) 

— oo 

The correlation function, which we shall study in this paper, is defined by 

tr (e''^ij{0,0)ij^{x,t)) 

(V'(0,0)V^(x,t))T = ^ ^ ^. (1.9) 

tr e T 

In the previous paper Q we obtained the determinant representation for this correlation func- 
tion. In this paper we shall derive completely integrable equations, starting from the determinant 
representation. The plan of this paper is the following. In Section 2 we shall remind the reader 
the determinant representation and definition of dual fields. In Section 3 we introduce new Hilbert 
space and rewrite the kernel of the integral operator in the canonical form. In Section 4 we define 
the resolvent of the integral operator. Section 5 is devoted to the construction of the Lax represen- 
tation. In Section 6 we find the logarithmic derivatives of the Fredholm determinant and obtain the 
completely integrable equation describing the correlation function. In Section 7 we summarize the 
main results. In Appendix A we find some identities for potentials. We give the treatment of the 
quantum nonlinear Schrodinger equation as a continuum generalization of the classical equation in 
Appendix B. Appendix C is devoted to the free fermion limit. 



2 Determinant representation for the correlation function 

Our starting point is the determinant representation for the temperature correlation function of local 
fields obtained in 
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(V(0,0)V't(x,i))T = e-^"*(0| (G{x,t) + ^) 
det ^i+V-^Y 



det (/ - ^Kt) 



|0) 



(2.1) 

a=0 



Let us explain our notations. 

We begin by the numerator in the r.h.s. of ( |2.1D , is the Fredholm determinant of the integral 
operator I + V — (here / is identical operator: I(A,^) = 5[X — /x)). This operator acts on the 
real axis. The left and right actions on some trial function / are given by 

CXD 

{l + V--^Y)of{^) = f{X)+ J (^v{X,f,)-^Y{X,f,)^ /(^)dM, 

— oo 

(2.2) 

oo 

/(A)o(/ + y-^y) =/(^)+ J /(A) (^F(A,/.)-£y(A,^)) dx. 

— oo 

Here and hereafter we denote by the symbol "o" the action of integral operators on functions. 

The kernels of operators V and Y can be written in terms of auxiliary quantum operators — 
dual fields, acting in an auxiliary Fock space. One can find the detailed definition and properties of 
dual fields in Section 5 and Appendix C of Q. Here we repeat them in brief. 

Consider an auxiliary Fock space having vacuum vector |0) and dual vector (0|. Three dual fields 
ipiX), (pDiiX) and (PA2W acting in this space are defined as 

0^2 (A) = QAiW+PDiW, 

cPmW = Qd^W+PaAX), (2.3) 
V'(A) = qi,{X) +p^{X). 

Here p{X) are annihilation parts of dual fields: p{X)\0) = 0; q{X) are creation parts of dual fields: 
(0|g(A) = 0. Thus, any dual field is the sum of annihilation and creation parts. Nonzero commutation 
relations are (see Q) 

[PAAX),q^{n)] = In h{fi,X), 
bD2(A),gv(Ai)] = ln^(A,/i), 

ic 

[PV'(A),9A2(/")] = ln/i(/x,A), where h{X, n) = — (2.4) 

A — fJj -\- tC 

[Pt/.(A),gDi(/i)] = ln^(A,^), 
bv(A)>9v(^)] = ln[/i(A,/i)/i(^, A)]. 
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Recall that c is the coupling constant in (1.2). It follows immediately from ( p.4| ) that the dual fields 
belong to an Abelian sub-algebra 



[V(A),V'(/i)] = [V'(A),<Aa(^)] = m\)Aa{l^)] = 0, 



(2.5) 



where a,b = A2,Di. This property, in fact, permits us to treat the dual fields as some c-number 
functions. 

Let us define function Z{X,fj.): 



The kernel V is equal to 



F(A,/x) 



47r2(A-/i) 



du 



+ 



<t>A2 («) 



Z{u,u) \u — X — iO u — X + iO u — fj, — iO u — /x + iO, 



Integral operator Y is one-dimensional projector 



where 
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Y{X,fi)=P{X)P{fi), 



+ 



27r j Z{u,u) \u — n — iO u — n + iO 

— oo ^ 

xe^(«)+-(«)e"5('/'M+-(^))Z(7/,/x). 



Here functions 9{X) and r(A) are equal to 



9{X) 



1 + exp 



e(A) 



(2.6) 



(2.7) 



(2. 



(2.9) 



(2.10) 



r(A) = itX^ — ixX. 



(2.11) 



The Fermi weight 9{X) defines the dependence of the correlation function on temperature T. The 
energy of one-particle excitation e(A) is given in ( |1.7| ). The function t(A) depends also on the 
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distance x and the time t. All other functions entering expressions for ^(A, fi) and do not 

depend on x and t. 

We would like to draw reader's attention to the fact that formulae ( p.TD and ( |2.9| ) slightly differ 
from formulae (6.25) and (6.26) of Q. It is explained in the Appendixes C and D of [Q] how one can 
reduce formulae (6.25) and (6.26) to formula ( p.7| ) and ( |2.9| ). Thus, the integral operator I+V — -^Y 
is explained. 

The operator I — also acts on the whole real axis. Its kernel is given by 

^) = {^TW^) v^\/^- ^^-^^^ 

Finally the function G{x, t) in the r.h.s. of (|2.1D is equal to 



G{x,t) = ^ I e'^^^'^+^^^Uv. (2.13) 

2lT J 



Thus, we have described the r.h.s. of ( p.ll ). The temperature correlation function of local fields 
is proportional to the vacuum expectation in the auxiliary Fock space of the Fredholm determinant 
of the integral operator. The auxiliary quantum operators — dual fields — enters the kernels V and Y. 
Due to the property ( ^.5]) the Fredholm determinant is well defined. Our aim now is a description 
of the correlation function in terms of solutions of classical completely integrable equations. 

3 Vectors and operators of new Hilbert space 

Introducing functions E±: 

^ , 1 Z(u,X) ( e~<^^2W e-'^^i(") \ r— 

^gVW+T(M)+|(</>Di(A)+0A2(A)-^(A)-T(A))^ (3,1) 

E^{\\u) = — Z(n,A)e^('^°iW+'^'*2W-^W-^W),/^. (3.2) 
27r ^ 

The functions E± depend also on the distance x, the time t, the temperature T and the chemical 

potential /i, but this dependence as a rule is suppressed in the notation. One can rewrite expressions 

(2.7) and ( |2.9| ) for V{\,^) and P()u) in terms of these functions 



V(\^) = I du{E+i\\u)EU\u)-E(\\u)E+U\u)), (3.3) 

X — fi 
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Y{X,^i) = P{X)P{^i) = J dudvE+{X\u)E+{n\v). (3.4) 

— oo 

Using obvious equalities 

oo 

d^E+{X\u) = -%E+{\\u)-t S dvg{u,v)E_{X\v), 



we arrive at 



(3.5) 



d,E_iX\u) = f^-(A|n), 



•\2 °° 

dtE+{X\u) = '-^E+{X\u)+ J dv{iXg{u,v)-d^g{u,v))E^{X\v) 



(3.6) 



dtE^{X\u) = -^^_(A|n). 



Here 



g{u,v) =6{u-v)e'^^'''^+^^''l (3.7) 

The relations ( |3.5D and ( |3.6[ ) are important for a description of the correlation function in terms 
of solutions of completely integrable equations. In order to do this it is necessary to investigate 
properties of the integral operator I + V. Later on we shall show how one can take into consideration 
the contribution of the projector Y. 

In order to derive differential equations for the correlation function it is convenient to treat 
functions E± as components of vectors of some Hilbert space Tl, for example, rigged L2(— oo, oo)<8'i?2 
Let us introduce the bra-vector (£'^(A)| and the ket-vector \E^{X)) belong to the Hilbert space 
TC. Both of these vectors have two discrete components (corresponding to the space R2), which 
we shall denote with indices 1 and 2. In turn any discrete component has continuous "index" 
(corresponding to the rigged L2(— 00,0x0) space) which we shall denote as "u" (or "u", "if" etc.). So 

( E^{^^\u) ^ 

{eHX)\ = {e^{X\u),eHX\u)); \E^{^^)) = 

The definition of the scalar product is standard : 

00 

{E'^{X)\E''{^,))= J dn[Ef{X\u)E^{fi\u)+E^{X\u)E^{fi\u)). (3.8) 

—00 

On the contrary the products of type \E^{fi)){E^{X) \ as usual are operators in Hilbert space TC. We 
shall consider such operators below. 
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Let us identify 

Ef{ii\u) = E+{fi\u), E^iX\u) = -E.{X\u), 

E^{ti\u) = E^if,\u), E^{X\u) = E+{X\u). 



(3.9) 



Then one can rewrite the kernel of the integral operator V as 

X — fi 



V{X,f,) (3.10) 

Due to (|3.9| ) we have 

(^^(A)I^^(A)) =0 (3.11) 

and hence the kernel V{X, fi) is not singular in the point A = /i. 

The representation ( 3.10| ) is the canonical form of the kernels of the completely integrable integral 



operators. In all examples related to correlation functions, the kernels of integral operators can be 
presented in the form ( 3.10| ). The different realizations of space TC correspond to the concrete corre- 



lation functions. For example, in the free fermion situation (the coupling constant c goes to infinity) 
TC = Rn, where n is the number of fields |^, |p, |^, [|l^]. For equal-time correlation functions of 
penetrable bosons the representations of type ( p. 10 ) were constructed with Ti. = ^2(0, 00) i?2n in 



, [11| (see also Section XIV of P]). In the present paper we shall follow the method developed 



in the papers enumerated above. 

Operators acting in the space 7i are defined in the standard way. They have discrete and 
continuous indices: A = Aji:{u,v), j,k = 1,2; —00 < u,v < 00. We shall denote these 
operators with the sign "hat" in order to distinguish them from integral operators which we have 
denoted with the sign "tilde" . Action on vectors is given by 

2 °° 

A\E^{fi)) = Y, I Ajk{u,v)Ej^{fi\v)dv, 

(3.12) 

2 

(ii;^(A)|i = ^ J Ef{X\u)Ajk{u,v)du. 
On the contrary the integral operators of type I + V appear to be scalars relatively to the space Ti, 
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for example 



^f(A|n)+ / V{X,fi)E^iii\u)dfi 



{i + v)o\E^{^,)) 



V 



J 



The product of operators of type A is 



i(i)i(2)=^ f Af{u,w)A 
/=i 

The trace of operator is defined as usual 

OO 

tr^= J du^Aii{u,u) + A22iu,u) 

— OO 

In particular for operators of type | ...)(... [ we have 

tr!...)(...| = (...|...). 



Using the operator notations, one can rewrite the relations ( |3.5D , ( |3.6| ) in the form of linear 
partial differential equations : 



5,|i?«(A)) = L(A)|i?«(A)), 9,(i?^(A)| = -(i?^(A)|L(A), 
dt\E^{X)) = M{X)\E^{X)), dt{E^{X)\ = -{E^{X)\M{X), 



where 



L(A) = Xa + [g,a], 
M(A) = -AL(A) + d^g, 
and the operators a and g are equal to 



<7 = -^cJ3<5(u- = j ^ ^ \6{u-v), 



-1 

9 = -cr+9{u,v) = I g{u,v) 





(3.13) 
(3.14) 

(3.15) 
(3.16) 

(3.17) 
(3.18) 



Later we shall use the relations ( |3.13| ) and ( |3.14| ) in order to construct the nontrivial Lax represen- 
tation. 



4 Vectors (F^(A)|, \F^{X)) and resolvent 

Let us introduce the vectors {F^{X)\ and |F^(A)) belong to the same space H 

( Ffifi\u) ^ 

defining them as solutions of the integral equations 

|F^(A))o(/ + y) = \E^if,)). 

More preciously these formulae mean 

oo 

fL[X)+ J y(A,^)i^^(/i)d/i = i^/(A), 



OO 

Ffif^)+ J Ff {X)V{X,fi)dX = Ef{fi). 



Define the resolvent of the operator I — V as 
Obviously 

[i - R) o {E\^)\ = {F\X)\, 

\E^{X))o(i-R) = \F^{^)). 
Let us find the kernel of the resolvent. One can rewrite ( [4.4[ ) as follows 

oo 

V{X,^)- j V{X,v)R{v,^j)du = R{X,n). 

— oo 

Multiplying both sides of the last equality by A — /i we get 

oo 

(A-/i)F(A,^) - / V[X,u)[X-u + v-fi)R{v,^ji)dv = {X-ix)R{X,^i), 
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or 



{E'^{\)\E''{^x))- J {E'^{\)\E''{v))R{v,ix)du = (i +V) o{v - ^,)R{v,ii), 



or due to (|4 



{E\\)\F\^^)) = (i + v) o {u-^i)R{v,^i). 
Making I — R act on this equation from the left, we get 

{F\\)\F^{^,)) = {\-^,)R{\^,). 

Therefore we have come to the following expression for the resolvent kernel: 

A — ^ 

It is worth mentioning that this method of calculation of resolvent is a direct generalization of the 
method described in the Section XIV. 1 of 0. 

In the next Section we shall need the operator B (potential), defined as 

Rf\\\/TpL, 



B = J |F^(A))(^^(A)|dA. (4.8) 

— oo 

Obviously 

oo 

S = |F«(A))o(/ + y)o(F^(^)|= J \E^{^^)){F\^,)\d^,. (4.9) 

— oo 

The components of this operator are 

oo 

Bjk{u,v)= J Fj^{X\u)E^{\\v)d\, (4.10) 

— oo 

so 

^^f B^,iu,v) Buiu,v) ^ ^^^^^^ 

\ B2l{u,v) B22{u,v) 

The operator C defined in the similar way is also useful 

oo 

Rf\\\/TpL, 



C = J A|F^(A))(^^(A)|dA. (4.12) 

— oo 

The components of the operator C are 

oo 

C,k{u,v)= f XFf{X\u)E^{X\v)dX, (4.13) 



3 

— oo 

SO 



C= I ^''^'"'''^ ^''^'''''^ \. (4.14) 

C2l{u,v) C22{u,v) 
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5 The Lax representation 



In this Section we construct the Lax representation having nontrivial compatibihty condition. Na- 
mely, we establish the following relations 

d,\F^iX))=C{X)\F^{X)), (5.1) 

9i|F^(A))=>i(A)|F^(A)), (5.2) 

which are analogous to the relations (|1|) and (|1|). We shall prove that one can obtain the 
operators C{\) and using formulae ( |3.15D and (|3.16| ) with replacement g hy g + B. 

Let us calculate the derivative of V{X,^) with respect to x using formulae ( 3.13 ) 

A — jjj 

= -{E^{X)\a\E^ii,)). (5.3) 

Thus, from (^^) we get 

oo 

9,|F^(A))o(/ + y) - J |F«(A))(i?^(A)|a|i?«(/i))dA = L(^)|i?«(/i)), (5.4) 

— oo 

or using the definition of B 

a,|F«(A)) o(^i + v)-B- a\E^{t,)) = L(^)li^«(^)). 

Making I — R on this equality the right, we get 

a,|F«(A)) - B ■ a\F^{X)) = (L(/i) - L{X))\E^{fi)) o (/ - r) + L{X)\F^{X)). 

In the r.h.s. we have 

im - L{X))\E^{f,)) o (l - r) = -a / |i?«(^))(F^(M)!F«(A))d/i 



— oo 



= -a-B\F^{X)). 

Therefore 

d,\F^{X)) = {LiX) + [B,a]) |F^(A)), 

or 



d,\F^iX))=C{X)\F^iX)), (5.5) 
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where 

£.{X) = X&+[b,a], (5.6) 

and 

b = B + g. (5.7) 
In the same way one can obtain the similar formula for (F^(A)| 

5.(F^(A)| = -(F^(A)|£(A). 

Now let us turn to the derivative of \F^{X)) with respect to t. As before, we start with differen- 
tiation of the kernel V{X, /x) using formulae ( |3.14[) 

' A — /i 

(ii;^(A)|((A - /i)2(A) - Km - LjXm^in)) 
X — fj, 

= {E^{X)\L{X)\E^iK)+MEHm\E''{K) 

= -d,{E^{X)\-\Em)+f^{EHm\Em)- (5.8) 

Thus, from (^) we get 

oo 

dt\F^{X)) o[i + v)- j \F^\X))d^{E\X)\ ■ \E^{f^))dX 

—oo 

oo 

+^ I \F^iX)){E^{X)\a\Em)dX = [-tiLit,) + d,g]\E^{fi)). (5.9) 

— oo 

Comparing ( ^ ) with ( |5.4| ), we see that 

oo 

[ \F^{X)){E^{X)\a\E^{f,))dX + ^,L{f,)\E^{f,)) 



= a,.!F^(A))o(^/ + yj 

Substituting this formula into (|5.9| ) we find 

at|F^(A)) o (i + v) + a,.|F^(A)) o (/ + . 
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= J |F^(A))9,(S^(A)|-|i?«(/.))dA + a,5|i^''(A^)). (5.10) 

—oo 

Acting on ( |5.10| ) by the resolvent from the right we have 

dtlF'^iX)) + 5,.|F^(A)) o (/ + f) . ^ o (/ - i?) 

oo \ 

J \F^i^^))^.{E'^{^^)\d^ • |F^(A)) + a.5|F^(A)). (5.11) 



The second term in the l.h.s. of ( 5.1l| ) is equal to 

a^|F^(A)) o(i + v) .fio(i-R 



= 5,|F^(A)) o (/ + . - A + A) o (/ - i?j 
Aa,|F^(A)) - a,|F«(A)) o (/ + y) o (F^(/i)|F«(A)) 

oo 

Aa.|F^(A)) - I / a.|F^(M)) • {E^{f^)\df^ 1 |F«(A)). 



Therefore we arrive at 



oo 



or 



Here 



dt\F^iX)) = -Xd^\F^{X))+\J ^,\F^{^,))■{E^{^^)\d^,j -iF^iX)) 

+ ^ \F^{^,))-dAE^if,)\d,)j ■\F^iX))+d,g\F^iX)) 

= -Aa,|F«(A))+a,(^ + 5)lF^(A)), (5.12) 

5i|F^(A)) = ^(A)|F«(A)). (5.13) 

Mix) = -XC{X) +d^b, (5.14) 



or, using (|5.6D we can rewrite (|5.14D as 

M{X) = -X'^a - X[b,a] + d^b. (5.15) 
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In the same way one can obtain the similar formula for {F^{X)\ 

dt{F^{X)\ = -{F^{X)\M{X). 
Thus, we have constructed the Lax representation. 

6 The logarithmic derivatives of the determinant and the differen- 
tial equations 

Remind the reader that the operator describing the correlation function contains the projector Y: 

oo 

Y{X,n) = P{X)P{^) = J dudvE+{X\v)E+{n\u). 

— oo 

In the operator notations it can be written in the form 

dudvE+{X\v)E+{n\u) = {E^{X)\a-\E^{fL)), (6.1) 



where the operator &- is equal to 

' 



£7- 



1 

(This is the particular case of operator — its discrete components are constant functions of continuous 
indices u and v. The action of this operator on vectors is still given by ( ^.121) ). 
Let us calculate the derivative of det (^I + V — with respect to a 



a«iogdet(/ + y-£y 



0=0 



= -i-TV ((i?^(A)|a_|ii;^(/z)) o (/ - i?)) = -i-Tr ((i?^(A)|a_|F^(A)) 

OO 

= -—tr (^Ba^^ = -— / Bi2{u,v) dudv. 

— oo 

Here we denote by the symbol "Tr" the trace of the integral kernels. Hence, 

CO 

dudv 



dadet{I + V- -^Y 



0=0 



det[i+v) J ^i?i2(n,^). (6.2) 



-oo 
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Due to the definitions (3.18|), (pj) and identity 



oo 

1 



G{x,t) = —— / dudvg{u,v), (6-3) 
2tt J 

— oo 

we come to the equahty 

oo oo 

G{x,t) — —— / Bi2{u,v) dudv = ——— / bi2{u,v) dudv. 

— oo — oo 

Hence the correlation function is equal to 

e-iht det(i+v) y 

(V^(0,0)V^t(^^^))^ = _ (0| V J / bu{n,v)dudv\0). (6.4) 

27T det(/- 2^Kt)_{^ 

The logarithmic derivatives of the determinant of the operator ^1 + with respect to x and t 
are equal 

log det (l + = Tr (^d^V o (^i - rJ 



dt log det (i + v)= Tr {dtV o {l - R 

Using (|5.3D we have 

Tr(a^yo (i-R)) = -TY[{E^{\)\a\E^{^j)) o [l - R 



dX{{E^{X)\a\F^{X))) = -tr(^S-a^. 



and hence 



log det (l + = - tr ' • (6-5) 
In the same method one can compute the time derivative 

= Tt [{E^{X)\{a{2X + ^-X) + [g, a])\E^{^^)) o (/ - r) 

OO 

dA(2A(^^(A)|<T|F^(A)) + {E^{X)\[g,a]\F^iX)) 



dXd^,^{E^{X)\a\E^{^^)){FH^^)\F^{\)) 

ti(2C ■ a - B"^ - a + B ■ [g,a]) = tr ( 2(0 + 6 • 5) -a-P -a 
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Therefore 

9tlogdet(/ + F) =11(2(0 + b- g) ■ a -P ■ a^. (6.6) 

The second logarithmic derivatives of the determinant can be expressed in terms of matrix b 
only. Indeed, using ( |A.5| ) and ( |A.7| ) we get 



log det (/ + y) = - tr Aft, a] • 6 • (tV (6.7) 



a^9tlogdet(l + y) = -tT[d^b-b-a -b-dj-aj. (6.8) 

Thus, we have expressed the logarithmic derivatives of the Fredholm determinant in terms of the 
traces of operators b, B and C. It is worth mentioning that the second logarithmic derivatives 
depend on the operator b only. 

Now let us turn back to the Lax representation. The relations ( ^.5| ) and ( |5.13| ) 

9.|F^(A))=£(A)|F«(A)), 

dt\F^{\)) = M{\)\F\\)) 
should be compatible. It means that 

dtC - d^M + [C, M] = 0. 
Substituting ( ^.6D and ( |5.15| ) into the last formula, we get 



[dtb.a] - d^d^b + 



[b,a],dxb 



0. (6.9) 



Remark. More accurately the equation ( p. 91) is valid if there exists a sequence of complex numbers 
A„ such that the corresponding sequence of vectors |F-^(A„)) generates a basis of the space H. 
However, one can check that this condition is not necessary. Indeed, using identities for potentials 



from Appendix A (in particular formulae ( A. 3 ) and (A^)) it is possible to prove the equality ( pS 
directly, without using the compatibility condition. 

Thus, we have arrived at the following results. The second logarithmic derivatives of the Fredholm 
determinant are presented in terms of the operator b. On the other hand this operator satisfies the 
partial differential equation ( |6.9D . In turn this equation appears to be the compatibility condition 
of the auxiliary linear problem ( ^.5]) , (^.13 ). 



17 



7 Main results in components of operators b and C 

In this Section we summarize the main results obtained in the previous Sections. Till now we did 
not use the symmetry of the kernel V{X,iJ,) = V{fx, A) (|3.3| ). Using this property and the definitions 
j!9| ) one can get additional identity 

Bu{u,v) = -B22iv,u). (7.1) 



Due to the definitions ( |5.7D and ( 3.1^ ) we have 



bab{u, v) = Bab{u, v), for all a, b exept the case a = 1, 6 = 2, 

6i2(n, v) = Buiu, v) - g{u, v) 
It follows from (^), ( |A^ ) and the definition that 



(7.2) 



dxBn{u,v) = i J dwbi2{u,w)b2i{w,v), (7.3) 

— oo 

oo 

dxB22{u,v) = -i dwb2i{u,w)bi2{w,v). (7.4) 



Using identities ( |A.5| ), ( A. 7 ) and expressions for the logarithmic derivatives of the determinant 
T5|)- (|6.8D we have 

oo 

9a; log det (^I + y j ~ ^ J d-ubii{u,u), (7.5) 

— oo 
oo 

9i log det + ~ *y ^'"'^^(~^2i(it, 'y)fi'(^^) ^t) 

— oo 
oo 

+ i J du{C22iu,u) - Cii{u,u))^ , (7.6) 

— oo 
oo 

dxdxlogdet(l + = - J dudvbi2{u,v)b2iiv,u), (7.7) 

— oo 
oo 

dtdx\ogdet(j + — J dudv{dxbi2iu,v) ■ b2i{v,u) 

— CO 

- dxb2i{u,v) ■bi2{v,u)). (7.8) 
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The diagonal part of the compatibihty condition ( |6.9[) is equal to zero due to (|7.3|) and ( |7.4[ ). 
For the off-diagonal part of the compatibility condition we find 

-idtbi2{u,v) = -dlbi2{u,v) 

oo 

+2 j dwidw2bi2iu,wi)b2i{wi,W2)bi2{w2,v), (7.9) 

— CO 

idtb2i{u,v) = -dlb2i{u,v) 

oo 

+2 / dwidw2b2i{u,wi)bi2{wi,W2)b2i{w2,v). (7.10) 



— oo 



The system of partial differential equations ( [7.S| ) and ( [7.10 ) is continuum generalization of the classi- 
cal nonlinear Schrodinger equation. The second logarithmic derivatives of the Fredholm determinant 
( |7.7| ) and (7.8) give us the densities of the first and the second integrals of motion of this equation. It 
means that the correlation function possess the properties of a r-function of the generalized nonlinear 
Schrodinger equation. 



Summary 

The main goal of this paper was a description of the quantum correlation function in terms of 
solutions of the classical completely integrable equation. We have constructed such differential 
equations — ( |7.9| ), ( [7.10 ), which, in fact, are continual generalization of nonlinear Schrodinger equa- 



tion. The correlation function of local fields plays role of the r-function of these equations. In 
particular the second logarithmic derivatives of the Fredholm determinant ( |7.7| ) and (7.8) give us 



the densities of the first and the second integrals of motion. In a forthcoming publication we shall 
formulate the Riemann-Hilbert problem for the differential equation obtained in this paper. This 
permit us to evaluate the long distance asymptotic. 



Acknowledgments 

We wish to thank Professor T. Inami for usefull discussions. This work is partly supported by 
the National Science Foundation (NSF) under Grant No. PHY-9605226, the Japan Society for the 
Promotion of Science, the Russian Foundation of Basic Research under Grant No. 96-01-00344 and 
INTAS-Ol-166-ext. 



19 



A Identities for potentials 

Here we establish some identities between matrix elements of potentials B (or b) and C. We also 
introduce the operator D: 

oo 

D= f X^\F^{X)){E^{X)\dX. (A.l) 



Let us calculate the derivative of the operator B with respect to x. We shall use formulae ( [3.13 ), 
and dUI). We have 



dj= J dx{£{X)\F^iX)){E\X)\-\F^{X)){E^{X)\L{X)) 

— oo 

dX(iXa+[b,a])\F^iX)){E^iX)\-\F^iX)){E^{X)\{Xa + [g,&]) 
= [a,C] + [b,a]-B-B-[g,a]. (A.2) 

Therefore 

d,B = [a,C] + [b,a]-B-B-[g,a]. (A.3) 
Taking the trace of this equality we have 



tr(^a,Bj=0. (A.4) 
Here we used definition b = B + g. Multiplying ( [A.3| ) by a and taking the trace we have 

iT(d^B-a]={.i([b,a]-b-a]. (A.5) 



Here we used the property 

g^ = {gaf = {), see O). 

Using the same method one can calculate the derivative of the operator B with respect to t. To do 
it we need formulae ( |3l^ ), ( p6D , (|5l^ ) and ( |5l^ ). We have 



dtB= j dX[M{X)\F''{X)){E\X)\-\F''{X)){E'^{X)\M{X)] 

— CO 

= -[a,b]-[b,a]-C + C ■ [g, a] + d,b ■ B - B ■ d^g. (A.6) 
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Again multiplying by a and taking the trace we have 

iAdtB ■ <T ) = tr f [ct, C\ - g ■ a - {a ,C\ - a - h ^ dj) ■ B ■ a - B ■ d.j,g ■ a 



(Here we used cyclic permutation under the sign "tr"). Substituting [o", C] from ( [A.3| ) we get after 
simple algebra 

tr(^ati?-(7^ =\.Y{dJ)-h-a -h-d.^h-'a^. {A.7) 

These identities were used for calculation of the second logarithmic derivatives of the Fredholm 
determinant in Section 6. One can also use the method described above in order to prove directly 
the equation (|6^) 



[dtb, a] - dxdxb + 



[b,a],dj) 



0. 



B Quantum equation as continual differential equation 

Consider the quantum nonlinear Schrodinger equation 

= -dli) + Icilj^ipip, (B.l) 

where ip{x, t) is annihilation operator in the bosonic Fock space. Let us take a matrix element of 
this equation between the ground state at finite density \Vt) and state with one hole (/io|- We shall 
have in mind that the momentum of the hole /cq is fixed. 

idt{ho\m = -dl{ho\m) + 2c(/io|V'^|/io, hi){hiMi'\h2){h2\m)- (B.2) 

Here /i2 is a hole with momentum k2 and hi is a hole with momentum ki. In the last term in the 
r.h.s. of ( |B.2| ) we have to integrate with respect to ki and /c2. This shows that quantum nonlinear 
Schrodinger equation is equivalent to the classical continual differential equation. The same type of 
equations describes quantum correlation functions (see (|7.9|) , (|7.10| )). 



C Free fermion limit 



In the free fermion limit c — > oo function E^{\\u) (3.2) does not depend on continuous variable u. 
Indeed, in this limit we have /i(A,^) = (A — /i + ic)/ic 1, hence all commutators ( |2.4| ) are equal 
to zero. Thus, one can put all dual fields equal to zero in expressions for functions E± 



E^{X\u) = ^{^—K + ^—^e<-^~'-^<'\ (C.l) 
2-K \u — A + zO u — X — J 

E_{\\u) = -e-5-(^)./^. (C.2) 
vr * 
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Hence matrix elements of the operator B possess the following properties: B2i{u,v) does not 
depend on u and v; Bii{u,v) depends only on the first argument u; B22{u,v) depends only on 
the second argument v; Bi2{u,v) depends on both arguments u and v. Analogous properties have 
matrix elements of the operator C. 

One can make the replacement 

B2i{u,v) ^ B21, C2iiu,v) ^ C21, 

00 00 

/ duBii{u,v) ^ Bii, J duCii{u,v) ^ Cii, 

—00 —00 

00 00 

/ dvB22iu,v) ^ B22, J dvC22{u,v) C22-, 

—00 —00 

00 00 

/ dudvBi2{u,v) B12, J dudvCi2{u,v) — > C12. 



Here functions Bab and Cab are scalar potentials. In terms of these potentials one can rewrite the 
correlation function, the logarithmic derivatives of determinant and the partial differential equations 
in the form 

—iht 

(V(0,0)V^(x,t))r = -^6i2det(/ + y). 



The logarithmic derivatives of the determinant 

d^\ogdei(l + V) = ibii, 

atlogdet(/ + y) = + C22 -&21G), 

9i5a;logdet(/ + = -612621, 

dtd^ log det(/ + = i{dxbi2 ■ 621 - <9a;62i • 612). 

The differential equations. 
(One should integrate ( [7.9^ with respect to u and v) 

-idtbu = -92612 + 26^2621, 

idtb2i = -9^621 + 26^1612. 

Here bab = Bab for all a and b except a = 1, 6 = 2; 612 = -B12 — G. These results reproduce the 
results of H up to notations and rescaling the distance x and the time t. 
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